Transient field fluctuations effects in d+Au and Au+Au collisions at 
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The effect of fluctuations of strong color electric fields (SCF) on the baryon production in d+Au 
and Au+Au collisions at 200A GeV is studied in the framework of the HIJING/BB v2.0 model. 
It is shown that the dynamics of the production process deviates considerably from calculations 
based on Schwinger-like estimates for homogeneous and constant color fields. An increase of the 
string tension from kq = 1 GeV/fm, to in medium mean values of 1.5 to 2.0 GeV/fm and 2.0 to 
3.0 GeV/fm for d+Au and Au+Au, respectively, results in a consistent description of the observed 
nuclear modification factors -RdAu and -RauAu (that relates (rf) Au+Au and p+p collisions), and 
point to the relevance of fiuctuations of transient color fields. The differences between nuclear 
modification factors T^AuAu and 7?cp (that relates central and peripheral collisions) are also discussed. 
The measurement of multi-strange (anti)hyperons (H, fl) yields would provide a crucial test of the 
importance of SCF fluctuations at RHIC energies. 

PACS numbers: 25.75.-q; 25.75.Dw; 24.10.Lx. 
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I. INTRODUCTION 

While the phase transition from hadronic degrees 
of freedom to partonic degrees of freedom (quarks 
and gluons) in ultra-relativistic nuclear collisions is a 
central focus of recent experiments at the Relativistic 
Heavy Ion Collider (RHIC), data on baryon and hy- 
peron production has revealed interesting and unex- 
pected features that may be of novel dynamical origin. 
As an example, a baryon/meson anomaly [H, |2jUI3| 
is observed as a large enhancement of the baryon to 
meson ratio and as a large difference of the nuclear 
modification factor (NMF) between total charged Q 
and neutral pious (tt*^) [6j at moderate (intermediate) 
transverse momenta {2 < pT < 6 GeV/c). 

In previous papers 0, [1] we studied the possible 
role of topological baryon junctions, 0,0, and the ef- 
fects of strong color field (SCF) [1] in nucleus-nucleus 
collisions. We have shown, in the framework of HI- 
JING/BB v2.0 model, that junction-antijunction {J J) 
loops with an enhanced intrinsic transverse momen- 
tum kx ~ 1 GeV/ c, a default string tension kq = 1 
GeV/fm, and a diquark suppression factor (7qq=0.07) 
provide a partial explanation of the baryon/meson 
anomaly [7|. That model provides an alternative 
dynamical explanation of the data to recombination 
models Within HIJING/BB v2.0 0,[1 one of the 
main assumptions is that strings could survive and 
fragment [ill, [Hi , and in particular populate the mid 
to low pt range. In contrast, in the recombination 
picture [lOj ] or in the hydrodynamical approach ^IS] 
all coherent strings are assumed to become rapidly 
incoherent resulting in a rapid thermalization. 

In nucleus-nucleus collisions the color charge exci- 
tations may be considerably greater than in nucleon- 
nucleon collisions due to the almost simultaneous in- 
teraction of several participating nucleons in a row 
[III . [31 and could be important even in the case of few 



binary collisions. Molecular dynamics models have 
been used to study the effects of color ropes as an 
effective description of the non-perturbative, soft glu- 
onic part of QCD [H. 

Recently, the effects of gluon field generated in the 
wake of hard processes and through primordial fluc- 
tuations of the color charges in the nuclei have been 
investigated [TgI [TtI. [T1| . The physical situation imme- 
diately after the collision bears close analogy to string 
models of high energy collisions. It was shown that the 
initial electric and magnetic fields produced in high 
energy hadronic collisions are longitudinal and leads 
to a novel string-like description of the collisions and 
a large topological charge density after the collisions 
[fTj . In our model, HIJING/BB v2.0, there are lon- 
gitudinal electric fields induced by a collision, which 
subsequently decay by quantum pair production. The 
effects of longitudinal magnetic fields are taken into 
account by a specific (like Mercedes star) topology 
for junction anti-junction J J loops |l9l|. Strangeness 
enhancement Jfepi, [2l|, [2^, |23|, \M, HfTstrong baryon 
transport [1, |26| |. and increase of intrinsic fc^ [3 are 
all expected consequences of SCF. These are modeled 
in our microscopic models as an increase of the effec- 
tive string tension that controls the quark-anti-quark 
{qq) and diquark - anti-diquark (qqqq) pair creation 
rates and the strangeness suppression factors [3| . 

In previous studies [?', 'sj we have focused on lon- 
gitudinal rapidity distributions and transverse mass 
(or momentum) spectra of hadrons and hyperons in 
Au+Au collisions at RHIC energies. However, the 
benchmark test of microscopic models are p+p and 
p{d)+A data. The results from p+p collisions are 
used as a base line reference to obtain nuclear mod- 
ification factor Raa- Therefore, we explore further 
dynamical effects associated with long range coherent 
fields (i.e strong color fields, SCF), including baryon 
junctions 9] and loops [1,[23, [lH, in the framework 
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of HIJING/BB v2.0 with emphasis on the novel 
baryon and hyperon observables measured at RHIC 
in p+p and d+Au coUisions. Using this model we an- 
alyze baryon/meson anomaly, the particle species de- 
pendence, centrality dependence in d+Au and AuH-Au 
collisions, as well as the differences seen in nuclear 
modification factors Raa (that relates Au-f-Au and 
p+p collisions) and Rqp (comparing central and pe- 
ripheral collisions). Comparison of NMF in different 
collision systems should provide information on the 
hadronization mechanisms. 

The paper is organized as follows. In Sec. II we re- 
view the HIJING/BB v2.0 (with SCF) model. For clar- 
ity, we repeat here some of the basics out-line that has 
been already presented in Refs. f?, "s*]. In Sec. Ill we 
discuss theoretical predictions in comparison with re- 
cent RHIC experimental data. Finally summary and 
conclusions are given in Sec. IV. 

II. HIJING/BB V2.0 MODEL 

A. Junction anti-Junction Loops 

HIJING is a model that provides a theoretical frame- 
work to extrapolate elementary proton-proton multi- 
particle phenomena to complex nuclear collisions as 
well as to explore possible new physics such as en- 
ergy loss and gluon shadowing [30|. Our analysis are 
performed in the framework of the HIJING/BB v2.0 
model that is based on HIJING/BB vl.lO [13, [Hi and 
HIJING [13]. 

In HIJING [lOl the soft beam jet fragmentation is 
modeled by diquark-quark strings as in [3l| with gluon 
kinks induced by soft gluon radiations. The mini-jet 
physics is computed via an eikonal multiple collision 
framework using pQCD PYTHIA 7.3 to compute the 
initial and final state radiation and hard scattering 
rates. In PYTHIA the cross section for hard parton 
scatterings is enhanced by a factor K = 2 va. order 
to simulate high order corrections. HIJING extends 
PYTHIA to include a number of new nuclear effects. 
Besides the Glauber nuclear eikonal extension, shad- 
owing of nuclear parton distributions is modeled. In 
addition dynamical energy loss of the (mini)jets is 
taken into account through an effective energy loss, 
dE/dx [30, 32]. 

In HIJING/BB vl.lO [l^l the baryon junction mecha- 
nism was introduced as an extension of HIJING/B [s^ 
in order to try to account for the observed longitudinal 
distributions of baryons(i3) and anti-baryons(i?) in 
proton nucleus {p+A) and nucleus- nucleus {A+A) col- 
lisions at the SPS energies. However, as implemented 
in HIJING/BB vl.lO the junction loops fails to account 
for the observed enhanced transverse slope of anti- 
baryons spectra at moderate pr in A+A [s^l- This 
is due to limitations of the pr algorithm adopted in 
version 1.10 that includes kinematic constraints that 
worked to oppose the predicted enhancement in the 
baryon junction loop In HIJING/BB v2.0 we re- 
placed that algorithm with one that implements J J 
loops pt enhancement directly. This is done by spec- 



ifying the intrinsic (anti)diquark pT kick in any stan- 
dard diquark-quark string (qq-g) that contains one 
or multiple J J loops. In addition, we introduced a 
new formula (see Eq.(l) below) for generating the 
probability that a given diquark or anti-diquark gets 
an "enhanced px kick" from the underlying junction 
mechanism. We emphasize that while there is very 
strong evidence from a variety of data that the source 
of the observed baryon pT enhancement arises more 
naturally from collective hydrodynamic flow [l3|, el- 
liptic flow of heavy hyperons (33 | argues strongly for a 
dominant partonic collective flow as the origin of the 
baryon anomaly. Most of the initial baryon radial px 
could theoretically arise from the production mecha- 
nism, while its elliptic deformation would arise from fi- 
nal state interactions. Nevertheless, our purpose here 
is to explore more fully without the kinematic limita- 
tions of version 1.10, the theoretical problem of how 
much of the baryon anomaly could be due to the pos- 
tulated baryon junction dynamics [^ at RHIC ener- 
gies and how well is described (p) c?-|-Au collisions with 
emphasis on strangeness production. The details of 
this new implementation of J J loops are described 
below. 

Multiple hard and soft interactions proceed as in HI- 
JING. Before fragmentation, we compute via JETSET 
[31] the probability that a junction loop occurs in the 
string. A picture of a juction loop is as follows: a color 
flux line splits at some intermediate point into two flux 
lines at one junction and then the flux lines fuse back 
into one at a second anti-junction somewhere further 
along the original flux line. The distance in rapidity 
between these points is chosen via a Regge distribu- 
tion as described below. For single inclusive baryon 
observables this distribution does not need to be spec- 
ified. 

The probability of such a loop is assumed to in- 
crease with the number of binary interactions, Uhits, 
that the incident baryon suffers in passing through the 
oncoming nucleus. This number depends on the rela- 
tive and absolute impact parameters and is computed 
in HIJING using the eikonal path through a diffuse 
nuclear density. 

We assume that out of the non single diffractive NN 
interaction cross section, (Tnsd=o'in — Csdf, a fraction of 
events fjj = o'jj/{aia — cTgdi) excite a junction loop 
(where CTsdf is single diffractive cross sections). The 
probability after rihits that the incident baryon has a 
J J loop is: 

PjJ = 1 - (1 - (1) 

We take a jj— 17 mb, tTsdf ~ 10 mb, and a-m ~ 42 mb 
for the total inelastic nucleon-nucleon cross section at 
nucleon-nucleon (NN) centre of mass (cm.) energy 
y/JNN— 200 GeV. These cross sections imply that at 
y/SNN=200 GeV, a junction loop occurs in p+p colli- 
sions with a rather high probability 17/32 « 0.5 and 
rapidly approaches 1 in collisions. In p-l-S" where 
Tihits ~ 2 there is an 80% probability that a junction 
loop occurs in this scheme. Thus the effects of loops is 
taken here to have a very rapid onset and essentially 
all participant baryons are excited with J J loops in 
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A+A collisions at RHIC. However, the actual proba- 
bility is modified by string fragmentation due to the 
threshold cutoff mass Mc = 6 GeV, which provides 
sufficient kinematical phase space for BB pair pro- 
duction. 

We investigated the sensitivity of the results to the 
value of parameters J J and found no significant vari- 
ation on pseudo-rapidity distributions of charged par- 
ticles produced in Au-|-Au collisions at ^/snn=200 
GeV, assuming cross sections ajj, 15 mb < ajj < 25 
mb and a cutoff mass Mc, 4 GeV < < 6 GeV. Light 
ion reactions such as, p+A or d-f-Au should show more 
sensitivity to ajj. In HIJING/BB v2.0, we introduced 
the possible topology with two junctions [l^ and a 
new algorithm where J J loops are modeled by an en- 
hancing diquark pT kick characterized by a gaussian 



width of CTq 



/ • Cqq, where / is a broadening factor. 



and (Tqq=0.350 GeV/c (consistent with PYTHIA |3l| 
default value). In Ref. [7|], we concluded that a value 
/ = 3 best reproduced the (anti)proton px spectra 
in Au-|-Au collisions at y/sNjf =200 GeV. This imple- 
mentation of the J J loops mechanism marks a radical 
departure from that implemented in HIJING/BB vl.lO. 
While the above procedure allows baryon anti-baryon 
pairs to acquire much higher transverse momentum in 
agreement with observation, the absolute production 
rate also depends on the diquark/quark suppression 
factor 7qq (see Sec. II B). The factor / modifying 
the default value a-qq=0.350 GeV/c may depend on 
beam energy, atomic mass number (A), and central- 
ity (impact parameter). However, we will show that a 
surprising good description of a variety of observables 
is obtained with a constant value, /=3. The sensitiv- 
ity of the theoretical predictions to this parameter is 
discussed in Sec. HI. Finally, we remark that baryon 
anti-baryon correlations studies in p+p and p{d)+ Am 
collisions at RHIC energies could eventually help us 
to obtain more precise values of J J loops parameters 
(mainly: (Tjj^ Regge intercept a{0), and parameter / 

)• 

Phcnomenological descriptions are currently based 
on Regge trajectory which gives the appropiate rela- 
tionship between the mass M of the hadrons and its 
spin Jg! Js = a{0) + AP, where a(0) is the Regge 
intercept, and is the Regge slope. The value of the 
Regge slope for baryons is as — 1 GeV~^ [1^ that 
yields a string tension (related to the Regge slope, 
kq = l/27ras) kq «1.0 GeV/fm. The multi-gluon 
exchange processes dominated by Pomeron exchange 
in high energy nucleus-nucleus collisions could be de- 
scribed by a Regge trajectory with a smaller slope 
as ' ~ 0.45 GeV^^ [35] , leading to an increase of string 
tension to k ~ 2 kq- 

The contribution to the double differential inclusive 
cross section for the inclusive production of a baryon 
and an anti-baryon in NN collisions due to J J ex- 
change is given by [1, [13] : 



Ci3se("(°)-i)l^«-^^l (2) 



■ (PpBfPpB 

where Cqq is an unknown function of the 



transverse momentum and Mq + P + B (junc- 
tion-t-Pomeron-fbaryon) couplings [13]. The pre- 
dicted rapidity correlation length (1— a(0))~-^ depends 
upon the value of the Regge intercept a(0). A value 
of a(0) ~ 0.5 ]9] leads to rapidity correlations on the 
scale \yB — Usl 2, while a value a(0) ~ 1.0 [sH] 
is associated with infinite range rapidity correlations. 
Thus, it is important to study rapidity correlations at 
RHIC energies where very high statistics data are now 
available. 



B. Strong Color Field within HIJING/BB v2.0 

In the case of quark-gluon plasma (QGP) creation 
it is necessary to modify the dynamics of particle vac- 
uum production at short time scales and the abun- 
dance of newly produced particle may deviate consid- 
erabl y fr om the values obtained assuming a constant 
field [33, HI] . Two possible processes that lead to an 
increase of strangeness production within the frame- 
work of the Schwinger mechanism are: i) an increase 
in the field strength by a modified string tension k 
[Hj [13, HI] , or ii) a drop in the quark masses due to 
chiral symmetry restoration [iO, [4l|, [l^, \^ ■ A spe- 
cific chiral symmetry restoration could be induced by 
a rapid deceleration of the colliding nuclei (i^ . 

For a uniform chromoclcctric flux tube with field 
(E) the probability to create a pair of quarks with 
mass m, effective charge ( Ceff) , and transverse momen- 
tum {pt) per unit time per unit volume is given by 
i45|]: 



P{pt) d^PT 



ecsE\ 

47^3 



ln< 1— cxp 



7r(m^ +Px) 



\ee«E\ 



(fpT 



(3) 

The integrated probability (Pm), when the leading 
term in Eq. 3 is taken into account, is given by: 



(eeff^ 

47r3 



1 

— exp 



TT m n 



eesE\) 



(4) 



where each term in the sum corresponds to produc- 
tion of n coherent pairs, and E is an homogeneous 
electric field, and K~\ecfiE\ is the so called string 
tension. We note, thal^ Pm reproduces the classi- 
cal Schwinger results [46i] , derived in spinor quan- 
tum electrodynamics (QED) for describing positron - 
electron (e+e~) production rate. A sizable rate for 
spontaneous pair production requires ^^strong electric 
fields", and Jeeff-El/m^ > 1. 

Recently, non-perturbative gluon and quark- 
antiquark pair production from a constant chromo- 
electric field with arbitrary color index via vacuum 
polarization have been also investigated in a QCD for- 
malism [ZTj . Although the pt dependence of the rate 
of production is different because of the presence of 
the nontrivial color generators in QCD, the integra- 
tion over Pt, reproduced also Schwinger's result for 
total production rate used here to estimate the sup- 
pression of heavier fiavors. 



4 



The real fields emerging in heavy-ion collisions, 
could be inhomogeneous, with a space-time depen- 
dence. Up to now, no reliable and universal method is 
available for calculating pair production rates in inho- 
mogeneous electric fields (for a review see Ref. (isj). 
Different theoretical methods, such as functional tech- 
niques [49! or kinetic equations 'st', ^3§\ have been de- 
veloped. It was shown [49] that the above formula is a 
very good approximation in the limit of sufficiently in- 
tense fields if the field does not vary appreciably over 
space-time distances of less than m~^{m'^ /\ecsE\). In 
the case of a QGP creation in heavy-ion collisions, 
the characteristic time of the field variation is esti- 
mated to be of order of 1 fm/c and we have to 
consider a time-dependent homogeneous field. Tak- 
ing a soliton-like field, the production probability for 
strange quarks is 1.5 times larger than in stationary 
case |37|, and therefore, the abundance of newly pro- 
duced multi-strange particles may considerably devi- 
ate from the values obtained for the constant field. 

In HIJING/BB v2.0, as in HIJING, no attempt has 
been made thusfar to study possible modifications of 
Lund string fragmentation by back reaction effects as 
discussed in reference [soj . These effects, important 
for consideration of final state interaction kinetic the- 
ory, go beyond the scope of the present phenomeno- 
logical investigation. 

A value for the strength of the external field, based 
on an estimate of initial energy density of 50 GeV/fm'^ 
at RHIC 'Ml, is ~ 3.16 GeV/fm. Due to time de- 
pendence and fluctuations of the chromolectric field at 
the initial stage of the collision we may even expect 
higher values. It has been suggested that the mag- 
nitude of a typical field strength at maximum RHIC 
energies might reach 5-12 GeV/fm [26| . 

In general in microscopic string models the heav- 
ier flavors (and diquarks) are suppressed according 
to Schwinger formula (for homogeneous strong color 
field, E) gl: 



7Q 



PjQQ) 
Pirn) 



exp 



Q 



ml) 



(5) 



where n = |eoffi?| is the string tension] mg is the 
effective quark mass; Q~s for strange quark; Q=qq 
for a diquark, and q~u, d are the light nonstrange 
quarks. 

The main parameters of QCD, the coupling 
strength q;qcd and the quark masses, need to be de- 
termined precisely. However, present estimates 
of the current quark masses range from: to„ =1.5-5 
MeV; TOd=3-9 MeV, and TOs=80-190 MeV. For di- 
quark we consider mqq=450 MeV [s^. Taking for 
constituent quark masses of light non-strange quark 
Mu,d= 230 MeV, strange quark A/,=350 MeV Q, 
and diquark mass Mqq=550 ± 50 MeV as in Ref. [53|, 
it is obvious that the masses of (di)quark and strange 
quark will be substantially reduced at the chiral phase 
transition. If the QGP is a chirally restored phase of 
strongly interacting matter, in this picture, the pro- 
duction of strange hadrons will be enhanced In 
this possible decrease of the strange quark mass 



would lead to a similar enhancement of the suppres- 
sion factors, obtained (in microscopic models) by an 
increase of string tension [s^, EH . If the strange quark 
mass is reduced from 1/^=350 MeV to the current 
quark mass of approximately rus ~ 150 MeV (the 
actual values are mg =80-170 MeV, see Ref. [52]), we 
obtain from Eq. 5 for strangeness suppression fac- 
tor 7^ « 0.70. For a string tension increase from 
Ko=1.0 GeV/fm to k=3.0 GeV/fm, we obtain an iden- 
tical value 7] « 0.69. Moreover, if we consider that 
Schwinger tunneling could explain the thermal char- 
acter of hadron spectra and that, due to SCF effects, 
the string tension value k fluctuates, we can define an 
apparent temperature as function of the average value 
of string tension (< n >), T = ^< k > /2tt [ssI. [56l|. 

There is a debate in the study of qqq system on the 
relative contribution of the A-like geometry and the 
F-like geometry fl9l, ^Et\ , and on the stability of these 
configurations for the color electric fields ^] . In both 
topologies we expect a higher string tension than in an 
ordinary qq string {ky = a/Sko and ka — (3/2) kq)- 
It was shown [58| that the total string tension has 
neither the Y nor the A-like value, but lies rather 
in-between the two pictures. However, the Y config- 
uration appears to be a better representation of the 
baryons. If two of these quarks stay close together, 
they behave as a diquark [57| . In a dual superconduc- 
tor models of color confinement for the y-geometry, 
the flux tubes converge flrst toward the centre of the 
triangle and there is also another component which 
run in the opposite direction (mercedes star). They 
attract each other and this lower the energy of the 
F-configuration (Tgj . 

Our calculations are based on the assumption that 
the effective enhanced string tension (k), in both, ba- 
sic ropes (g" — (p) and junction ropes (g" — q" — q") 
are the same. For elementary n strings and junctions 
this ansatz is supported by baryon studies [59|] ■ A dif- 
ferent approach to baryon production without baryon 
junctions has been proposed in ^23| where SCF from 
the string fusion process can lead to {qq)(, — {qq)e with 
about a doubling of the string tension. Both types of 
SCF configurations may arise but predict different ra- 
pidity dependence of the valence baryons. 

Following the discussions above, we take into ac- 
count SCF in our model by an in medium effective 
string tension k > kq, which lead to new values for 
the suppression factors, as well as the new effective in- 
trinsic transverse momentum /cy . This includes: i) the 
ratio of production rates of diquark to quark pairs (di- 
quark suppression factor), 7qq = P(qqqq)/P(gg), ii) 
the ratio of production rates of strange to nonstrange 
quark pairs (strangeness suppression factor), 7^ = 
P{ss)/ P{qq), iii) the extra suppression associated 
with a diquark containing a strange quark compared 
to the normal suppression of strange quark (7^), 7us = 
(P(usus)/P(udud))/(7s), iv) the suppression of spin 1 
diquarks relative to spin ones (appart from the fac- 
tor of 3 enhancement of the former based on counting 
the number of spin states), 710, and v) the (anti)quark 
(cr^' = y^n/KQ ■ Gq) and (anti)diquark (cr^'q = ^/k/kq- 
f ■ (Tqq) gaussian width. These parameters corre- 



5 



spond to 7qq=PARJ(l), 7^=PARJ(2), 7us=PARJ(3), 
7io=PARJ(4), and tJqq = aq= PARJ(21) of the JET- 
SET subroutines [sij. The values of the main param- 
eters sets used in our analysis (i.e /, cr,, as well as the 



suppression factors 7qq, 
the Appendix. 



7s, 7us, and 710) are found in 



III. NUMERICAL RESULTS 

A. Transverse momentum spectrum 

1. Identified particles spectra from p+p collisions at 
y/sWN=200 GeV 

We will concentrate our discussions on species de- 
pendence of the nuclear modification factors (NMFs) 
i?AA and i?cp in li+Au and Au-|-Au collisions at 
^/sNN=2{)Q GeV. Raa is the ratio of the heavy-ion 
yield to the pp cross section normalized by the mean 
number of binary collisions (< A^bin >), while Rev 
is the ratio of scaled central to peripheral particle 
yield. They are defined as in Ref. By compar- 
ing the yields in p{d)+A and A+A collisions to that 
from p+p collisions, with a scaling factor to take into 
account the nuclear geometry, one can test the as- 
sumption that nucleus-nucleus collision is a simple 
superposition of incoherent nuclcon-nucleon scatter- 
ing and explore possible nuclear effects (e.g., shad- 
owing, quenching, SCF). The relevance of baseline 
p+p hard pT spectra for understanding high-energy 
nucleus-nucleus physics is discussed in Ref. [60|. The 
precision in the first results on the NMFs i?AA, were 
limited by the uncertainty in the parameterization of 
the p+p reference spectrum used in obtaining i?AA- 
This limitation was partly overcome when new data 
on p+p, d+Au and Au+Au collisions at ^/sJfN=200 
GeV using the same experimental setup were obtained 
[6l|,|62, 63, 64, 65, 66, 67, 68, 69, .7fl,[lll[7|. 

Before discussing our results on the NMFs that 
are scaled yield ratios, it is important to show how 
well the model describe p+p collision data. Fig. [T] 
presents a comparison of experimental transverse mo- 
mentum spectra |6^, H^l of positive pious (7r+), kaons 
protons (p), (multi)strange particles and their 
anti-particles, with the predictions of regular HIJING 
(dotted histograms) and HIJING/BB v2.0 (solid his- 
tograms). HIJING results describe rather well the 
data for produced "bulk particles" (i.e. pious, pro- 
tons, and kaons). We would have expected that 
other (multi)strange particles to also be well described 
by the standard fragmentation included in HIJING, 
since the universality of the fragmentation process be- 
tween positron-electron {e^ + e~) and p+p collisions 
has been confirmed [t^. However, regular HIJING 
underestimate significantly the spectra of lambda and 
cascade particles (Fig. [T]c,d). 

HIJING/BB v2.0 results include the contribution of 
a broadening of the intrinsic momentum (fc^), simu- 
lated by an increase of the string tension from the de- 
fault value Ko=l-0 GeV/fm (Tab. H Set 2) to ki=1.5 
GeV/fm (Tab. [H Set 3). Such parametrization is 
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FIG. 1: (Color online) Comparison of HIJING/BB v2.0 
predictions for the invariant yields of pions, kaons (scaled 
by 10~^), protons (scaled by 10~^) and their anti-particles 
(upper panels) and (multi)strange and their anti-particles 
(lower panels) as function of px from non-single diffractive 
(NSD) p+p coUlisions at ^5^=200 GeV. The rapidity 
range was —0.5 < y < 0.5. Results with ki=1.5 GeV/fm 
(solid histograms; Tab. I, Set 3) and with ko=1-0 GeV/fm 
(dotted histograms; Tab. I, Set 1) are shown. The data 
are from STAR [to . [13] • Only statistical error bars are 
shown. 



supported by earlier experimental measurements [73| , 
showing that at center of mass energy ^/snn > 
100 GeV, collisions between protons and anti-protons, 
largely consist of more than a single parton-parton in- 
teraction. It is also supported by a recent PHENIX 
analysis of the properties of jets produced in p+p 
colhsions at ^/snn —200 GeV [ill, giving a value of 
< kr >=1.34 ± 0.04 ± 0.29 GeV/c, much larger than 
the naive expectation for the pure intrinsic parton 
transverse momentum based on nucleon constituent 
quark mass (i.e., < kx > ^ 0.300 GeV/c) [zi]. This 
new parameterization results in an increase of the 
yield of (multi)strange particles without affecting sig- 
nificantly results for the "bulk particles" and thus 
gives a simultaneous description of the non-strange 
and multi-strange sector. 

A related approach has been recently presented by 
STAR collaboration [6^. They introduced an im- 
proved description for strange particles using PYTHIA 
by increasing the K factor which quantifies the contri- 
butions of next-to-leading order (NLO) effects, from 
K=2 to K=3. We test this suggestion using our model 
and observes that such increase of i^T results in an over- 
estimation of pions, kaons, and protons yields. In our 
calculations we keep K to its default value in PYTHIA 
and HIJING i.e., K^2 for all collisions. 



2. Identified particles spectra from d+An collisions at 
y/sWN=200 GeV 

Initial-state nuclear effects are present in both 
d+Au and Au+Au collisions, while final state effects 
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are expected to contribute only in Au+Au collisions. 
Thus, effects from the initial state are best studied 
through a "control" experiment such as rf+Au colli- 
sions [TTI . [tsI . [tqI . Isoj . Multiple soft scattering of pro- 
jectile partons as they traverse a target nucleus may 
increase their transverse momentum before they un- 
dergo the hard scattering or subsequent to it, leading 
to an enhancement of the yield at moderate and high 
Pt compared to p+p collisions - called the Cronin ef- 
fect [8l|, [82]. This enhancement is expected to have 
some particle mass dependence and to be stronger 
for heavier particles |82l |. The Cronin enhancement 
was addressed in Refs. [H, [13, [H, HI]. However, 
these pQCD calculations can not predict the parti- 
cle species dependence observed in the data, as initial 
state parton scattering precedes fragmentation into 
the different hadronic species. An alternative expla- 
nation considering final state interactions is discussed 
in Ref. 

HIJING type models incorporate in addition to the 
soft and hard dynamics (discussed in Sec. II A) a 
simulation of soft multiple initial state collision effects 
[9l| . Excited strings are assumed to pick up random 
transfer momentum kicks in each inelastic scattering 
according to the distribution: 

g{k,) cx {{kl+pl){kl+pl){l + exp[{k,-p2)lp3])y^ 

(6) 

where kc represents are ntrinsic transverse momen- 
tum of the colliding partons. The parameters pi=0.1, 
P2 = l-4, P3=0.4 GeV/c are chosen to fit the low en- 
ergy multiparticle production [s^] . A flag in the code 
makes it possible to compute spectra with and with- 
out this effect. At RHIC energies the observed larger 
enhancement for protons and anti-protons in compar- 
ison with pions |72l ] requires new processes beyond 
the initial state multiple scattering. In addition to 
"Cronin enhancement" other known initial state ef- 
fects that could contribute include nuclear shadowing 
[lO,!^,!!], gluon saturation ^ga], J J loops, and SCF 
effects Ref. Q. 

In order to further understand the mechanisms re- 
sponsible for the particle pT spectra, and to sepa- 
rate the effects of initial and final partonic rescatter- 
ings we study the transverse momentum spectra at 
mid rapidity of identified particles in central (0-20%) 
(Fig. HI) and peripheral (40-100%) (Fig. O d+kn col- 
lisions. The data are compared to the predictions 
of our model assuming two values of the string ten- 
sion, K2=2.0 GeV/fm (solid histograms) and ki— 1.5 
GeV/fm (dotted histograms). All calculations in- 
clude nuclear shadowing but do not take into consid- 
eration "jet quenching", since d-|-Au collisions create 
a "cold nuclear medium". The agreement above 1 
GeV/c, is very good for peripheral (< iVbin >=4.5) 
and minimum bias (< A^bin >=7.8) events (not shown 
here). Within our model the observed "Cronin en- 
hancement" for identified particle spectra at moder- 
ate p-r, as well as species dependence, are satisfactorily 
predicted at mid-rapidity due to a subtle interplay of 
initial state effects. 




4 2 4 

p, (GeV/c) p, (GeV/c) 



FIG. 2: (Color online) Upper part: comparison of HI- 
JING/BB v2.0 predictions for the invariant yields of pi- 
ons, kaons (scaled by 10~^), protons (scaled by 10"'^) (left 
panel) and their anti-particles (right panel) as function 
of Pt for central (0-20%) rf+Au colUisions at ^sjvjv=200 
GeV. Lower part: spectra for strange (A) and multi- 
strange (H~, f2~) particles (left panel), and their anti- 
particles (right panel). The rapidity range was —0.5 < y < 
0.5. Results with K2=2.0 GeV/fm (solid histograms; Tab. 
II, Set 4) and with ki=1.5 GeV/fm (dotted histograms; 
Tab. II, Set 3) are shown. The data are from STAR [H]. 
Only statistical error bars are shown. 



However, the calculations somewhat underestimate 
the pion, proton and kaon yield in the case of cen- 
tral d+Au collisions (< A'bin > = 14.8). To quantify 
the competing mechanisms contributing to the par- 
ticle yield, we study in the above results, the con- 
tribution from "soft Cronin" and shadowing effects. 
Including "soft Cronin" through Eq. 6 results in an 
increase of the yield which is higher for protons (« 
40% for pT=2-3 GeV/c) than for pions (« 20% for 
Pt=2-3 GeV/c). On the other hand shadowing ef- 
fects, result in a roughly equal decrease of the yield of 
all species, i.e., « 25% in the same pt region. There- 
fore, these two effects act in opposite directions and 
partly cancel each other at low and moderate pr {< 5 
GeV/ c) . In addition, we can not improve the overall 
description of pt spectra, by modifying the parame- 
ters from Eq. 6, without destroying the satisfactory 
results at low (0.0-1.5 GeV/c). A further increase 
of mean values of string tension to values greater than 
2.0 GeV/fm, results in a stronger underestimation of 
the pions yields. The description of pt spectra for 
central collisions point to a mechanism not included 
in our model. 

The disagreement may be explained by considering 
hadronic rescattering in gold nucleus as in Multi Phase 
Transport (AMPT) model [o^] , or inelastic and elastic 
parton ladder splitting as in EPOS model |93|] . We do 
not include here a discussion of the particle production 
at forward rapidity, where possible gluon saturation 
effects (color glass condensate CGC) [s^ or parton 
ladder splitting [§3], may have to be considered. 
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FIG. 3: (Color online) Results for peripheral (40-100%) 
d+Au coUlisions at -^sjvjv=200 GeV. The histograms have 
the same meaning as in Fig. [5] The data are from STAR 
[g^ . Only statistical error bars are shown. 



In Section III B we will discuss the nuclear modifi- 
cation factor for (multi)strange particles. Therefore, 
predictions for A, S~, Q~ and their antiparticles are 
also shown in Fig. [2] and Fig. [3] (lower panels) for the 
same centrality selections. We note that in rf+Au col- 
lisions the sensitivity to in medium string tension de- 
pend on mass and strangeness content, resulting in 
a predicted significant increase of the yield of multi- 
strange hyperons (S and fl). 



3. Identified particles spectra from Au+Au collisions at 
^/s^= 200 GeV 

Figure [4] presents a comparison of the transverse 
momentum spectra of identified 7r+, and p and 
their anti-particles with the predictions of HIJING/BB 
v2.0 (upper panels) for central (0-5%) Au+Au colli- 
sions. These data are taken from Ref. [6§|. A similar 
comparison (lower panels) is also presented for As (for 
central 0-5%), S~s, ft^s (for central 0-10%) and their 
anti-particles. These data are taken from Ref. |71|. 
The shape of the measured spectra show a clear mass 
dependence. The protons and anti-protons as well as 
hyperons (S, fl) spectra have a pronounced shoulder- 
arm shape at low pr characteristic of radial flow. The 
results of the model are shown by the dotted (^2=2.0 
GeV/fm) and solid (^3=3.0 GeV/fm) histograms. In- 
troducing the new J J loops algorithm and SCF effects 
in HIJING/BB v2.0 results in a significant improve- 
ment (relative to the predictions without SCF [71) in 
the description of the pion, kaon and proton spectra 
at intermediate px for both values of k. Only a qual- 
itative description is obtained at low pT due to the 
presence of elliptic (hydro) and radial expansion [isj , 
not included in our model. 

The yields and transverse momentum slopes of 
(multi)strange particles at moderate pr are still un- 
derestimated for the lower value of string tension 



FIG. 4: (Color online) Upper part: comparison of HI- 
JING/BB v2.0 predictions for the invariant yields of pi- 
ons, kaons (scaled by 10~^), protons (scaled by 10"'^) (left 
panel) and their anti-particles (right panel) as function of 
Pt for central (0-5%) Au-|-Au coUlisions at y^sjvjv=200 
GeV. Lower part: spectra for strange (A) and multi- 
strange (H~, n~) particles (left panel), and their anti- 
particles (right panel). The rapidity range was —0.5 < 
y < 0.5. Results with K3=3.0 GeV/fm (solid histograms; 
Tab. II, Set 5) and with K2=2.0 GeV/fm (dotted his- 
tograms; Tab. II, Set 4) are shown. The data are from 
PHENIX ^ and STAR [73, [zl- Only statistical error 
bars are shown. 



^2=2.0 GeV/fm. Multi-strange hyperons yields and 
spectra, seems to favour the larger value K3=3.0 
GeV/fm. This suggest that multi-strange particles 
are produced early in the collisions, when tempera- 
ture is higher than those which characterize the "bulk 
particles" (i.e., tt, K, p). This points towards a dy- 
namical origin and could be explained as an effect of 
fluctuations of the transient strong color field at early 
time. Recently, it was also shown that "bulk parti- 
cles" seems to have a different temperature at kinetic 
freeze-out than hyperons (S, ^2) suggesting also that 
multi-strange baryons do not take part in "i/ie same 
collectivity" as tt, K, and p during the collision [o^loH- 



B. Nuclear Modification Factors 

To evaluate if this version of HIJING/BB describes 
the produced entropy, the predicted transverse mo- 
mentum spectra (left panels) and RdAu (right panels) 
of the total hadron yield for central (0-20%) d+Au col- 
lisions, are compared to data from Refs. [H, [t^] in 
Fig. The data could not be well described by as- 
suming only a broadening of the intrinsic kT from its 
standard value crqq=0.360 GeV/c (dotted histograms) 
to CTqq=1.08 GeV/c (i.e., including J J loops, dashed 
histograms). The introduction of SCF has a slight 
effect on the predicted nuclear modification factors 
RdAu, but results in a somewhat better agreement 
with data (solid histograms). Similar results are ob- 
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tained for minimum bias data. The data indicate at 
most a small variation with centrality of the factor / 
consistent with the broadening originating at the par- 
ton level. 
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FIG. 5: (Color online) Left panel: HIJING/BB v2.0 pre- 
dictions with SCF for pr spectra at mid-rapidity of total 
inclusive charged hadrons, for central (0-20%) rf-|-Au and 
p+p collisions. The calculations assume ^2=2.0 GeV/fm 
for d-|-Au (solid histograms; Tab. II, Set 3) and ni — 1.5 
GeV/fm for p+p (dot-dashed histograms; Tab. I, Set 3). 
Right panel: Nuclear modification factor RdAu- The re- 
sults obtained without J J loops and SCF (dotted; Tab. 
II, Set 1), including only J J loops (dashed; Tab. II, Set 
2), and including both effects J J loops and SCF (solid 
histograms. Tab. II, Set 3) are shown. The data are from 
STAR [6^, [t^. Only statistical error bars are shown. 

A slight discrepancy is seen in the description of 
observed px spectra at midrapidity for total inclusive 
charged hadrons yields in central (0-20 %) d+Au colli- 
sions. In order to understand this, we analyse charged 
hadron asymmetries (i.e., the ratio of the hadron yield 
in backward rapidity to forward rapidity intervals) as 
proposed in reference [11] . Since this observable is de- 
fined as a ratio, it is important to check the results 
for both numerator and denominator. Fig. [5] presents 
the yields in different pseudo-rapidity interval, back- 
ward (part a) and forward (part b), obtained within 
HIJING/BB v2.0 and HIJING models. The absolute 
Pt distribution are better described by HIJING/BB in 
comparison to the results of regular HIJING. 
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FIG. 6: (Color online) HIJING/BB v2.0 (with SCF) and 
HIJING vl.382 predictions for pr spectra of total inclu- 
sive charged hadrons, for central (0-20%) d+Au collisions, 
in backward ( —1.0 < rj < —0.5, left panel) and forward 
(0.5 < 7] < 1.0, right panel) pseudo-rapidity region. The 
calculations assume «:2=2.0 GeV/fm for ci-|-Au (solid his- 
tograms; Tab. II, Set 4) and ki=1.5 GeV/fm for dotted 
histograms; Tab. II, Set 3). The dashed histograms are 
the results obtained within HIJING vl.382, with default 
parameters. The calculations in both models are obtained 
with shadowing and without jet quenching. The data are 
from STAR IQo]. Only statistical error bars are shown. 

In Figure[71 we analyse the observed charged hadron 



asymmetries within both models for central (0-20%) 
and minimum bias d-l-Au collisions. HIJING/BB and 
HIJING models predict similar asymmetries. The cal- 
culated asymmetries are greater than 1.0, in quali- 
tative agreement with the observations and contrary 
to pQCD predictions 97]. Within our model, the 
results on asymmetry are also little sensitive to the 
strength of the assumed longitudinal chromoelectric 
field. For central (0-20%) rf-|-Au collisions, increas- 
ing string tension results in a decrease of the abso- 
lute yield (see Fig. 0, without changing the asym- 
metry significantly. There are still some differences 
between the pr dependence of the predicted and mea- 
sured asymmetry. A relatively good agreement is ob- 
tained for minimum bias d+Au collisions in the region 
(0.0 < \ri\ < 0.5), not shown here. This observable is 
very delicate and great care should be taken on its 
interpretation. 
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FIG. 7: (Color online) The ratio of charged hadron spectra 
in the backward rapidity to forward rapidity region for 
minimum bias (upper panel) and central (0-20%) (lower 
panel) d+Au collisions. The curves have the same meaning 
as in Fig. [6] The data are from STAR ^96]. Only statistical 
error bars are shown. 

Our main conclusion is that the small asymmetry is 
not well reproduced by either models. These observed 
differences may be related also to the disagreement 
in the calculated pr spectra of identified particles for 
central d+Au collisions discussed above (Sec. IIIA). 

Higher sensitivity to the new dynamics imple- 
mented in the HIJING/BB v2.0 is obtained by study of 
nuclear modification factors of identified particles. To 
better quantify possible effects of strong color field on 
particle production, results of RdAuiPr) for central 
(< iVbin >=14.8) d+Au coUisons where higher sen- 
sitivity to SCF is expected, are presented in Fig. [51 
Because of their dominance, the production of pions 
is only moderately modified when we assume an in- 
crease of the string tension value since the total en- 
ergy is conserved. Taking into account SCF effects 
(solid histograms) results in changes at moderate pt 
of w 15-20% of the predicted pion yield (Fig.[HJi). The 
pions yield in central d+Au collisions is enhanced rel- 



9 



ative to p+p collisions (i.e., R^j^^iPT) > 1) an effect 
reproduced by our calculation. The scaling behavior 
of sum of protons and anti-protons {p+p) is somewhat 
different than that of the pions. The (anti)protons are 
much more sensitive to J J loops and SCF effects, and 
as a result R^dAu \pT) is predicted to be higher than 
that of pions at moderate pt- This result is consistent 
with the data of Ref. [g^ , where this behavior was ex- 
plained by different Cronin enhancement for pions and 
protons. However, within our model it is not only the 
"Cronin effect" which modified Rd Am but an interplay 
of initial state effects with an important contribution 
coming from J J loops and SCF. In contrast to pions, 
(multi)strange particles show a high sensitivity to the 
presence of SCF. Kaons (Fig.[5)3) and lambda particles 
(Fig. [8ji), show an increase at moderate pt by a factor 
of approximately 1.5 and 3.0 respectively, relative to 
the calculation that do not include J J loops and SCF 
effects (dotted histograms) . As opposed to pions such 
an increase results in a predicted strong enhancement 
of the lambda yield relative to scaled binary collisions. 
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FIG. 8: (Color online) HIJING/BB v2.0 predictions for 
species dependence of nuclear modification factor {Rd Au) 
in central (0-20%) d-|-Au collisions at i/sjvjv=200 GeV. 
Shown are results for; (a) charged pions, (b) kaons, (c) 
inclusive p+p, and (d) inclusive A -I- A. The histograms 
have the same meaning as in right pan el of Fig.O The data 
are from STAR (filled symbols) H and from PHENIX 
(open symbols) [72]. Only statistical error bars are shown. 



Figure M presents the values of RdAu for peripheral 
(40-100%) collisions, where the mean number of col- 
lisions is very small, < A^bin >=4.5. Our calculations 
shows that SCF effects play an important role already 
when there is only a few number of binary collisions. 
The calculation including SCF describe well the ex- 
perimental results. 

The sensitivity to transient field fluctuations is 
shown in Fig. [TO] which present the predicted 
RdAuiPT) assuming K2=2.0 GeV/fm and ki—1.5 
GeV/fm for central (0-20%) d+Au collisions. For 
k;i=1.5 GeV/fm the calculation predicts a value of 
RdAu close to one and that decrease somewhat with 



the mass and strangeness content. Increasing k to 
K2—2.0 GeV/fm results in a larger enhancement rel- 
ative to p+p (i.e., RdAu > 1) and a reversing of the 
order of the nuclear modification factor as a function 
of the mass and strangeness content. In particular, the 
model predicts an enhancement in the multi-strange 
(anti)hyperon production at moderate p^, by a factor 
of roughly 2 relative to binary scaling for Ss and by a 
factor of approximately 3 for fls. 
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FIG. 9: (Color onhne) RdAu for peripheral (40-100%) 
d+Au collisions at y^sjvjv=200 GeV. The histograms have 
the same meaning as in right panel of Fig. \5\ The data are 
from STARffiUed symbols) ^ and from PHENIX (open 
symbols) [T^l- Only statistical error bars are shown. 



Preliminary data on hyperons A and S [2ll . |99|, 
seems to favour ki — 1.5 GeV/fm. However, due to 
large statistics and systematics error, no clear con- 
clusion could be drawn. Measurements of H. particles 
could help us to draw a final conclusion on the impor- 
tance of SCF fluctuations at top RHIC energy and to 
determine the effective value (k) of string tension, in 
d-l-Au collisions. 

We show in Sec. IIIA, that a new parametriza- 
tion of the p+p interaction is necessary to describe 
the identifled particle spectra especially for strange 
and multi-strange particles. This parameterization 
is different than that of PYTHIA mainly in assum- 
ing different constituent and current quark masses for 
diquark. Therefore the strength of SCF (and k) in- 
troduced in our previous work Q to best described 
-RauAu has to be modified because of this new p+p 
baseline. Figure [TT] show the sensitivity of NMF 
i^AuAu to transient color field fiuctuations for the p+p 
yield and for (multi)strange particles. The results as- 
suming K=K2=2.0 GeV/fm and k—K3=3.0 GeV/fm 
are shown for both central (0-5%) and peripheral 
(60-90%) Au-f Au collisions. In agreement with the 
data, the calculated p+p yield scale approximately 
with the number of binary collisions. In contrast, the 
observed i?AuAu for (multi)strange particles show an 
enhancement over p+p, which increase with the mass 
and strangeness content. A good descriptions of NMF 
i^AuAu for (multi)strange at moderate pT, is obtained 
only for the largest value of k, i.e., K3 (Fig. [TTk). 
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FIG. 10: (Color online) HlJING/BB v2.0 predictions in- 
cluding SCF with K2=2.0 GeV/fm (Tab. II, Set 4) (up- 
per panels) and with ki = 1.5 GeV/fm (Tab. II, Set 3) 
(lower panels) for central (0-20%) (left panels) and min- 
imum bias (right panels) li+Au collisions at -^/sjvjv=200 
GeV. Shown are results for: inclusive p+p (solid), inclu- 
sive A + A (dashed), H+ +H~ (dash-dotted) and Q,'^ + 
(dotted histograms). The data are from STAR [6^. Only 
statistical error bars are shown. 
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FIG. 11: (Color online) HIJING/BB v2.0 predictions in- 
cluding SCF with At3=3.0 GeV/fm (Tab. II, Set 5) (upper 
panels), and with K2=2.0 GeV/fm (Tab. II, Set 4) (lower 
panels), for central (0-5%) (left panels) and peripheral 
(60-90%) (right panels) Au+Au collisions at ^sjvjv=200 
GeV. The histograms have the same meaning as in Fig. [TO] 
The data are from STAR [11,163,111. Only statistical error 
bars are shown. 



A clear ordering with strangeness content, seen in 
the data, is not reproduced by lower value k=K2—2.0 
GeV/fm (Fig. [Tib). This clearly illustrates that in 
Au+Au collisions, fluctuations of the chromoelectric 
field are higher than in d+Au collisions, and affect es- 
pecially (multi)strange particle yields. The p, 7r+, A'+, 
and their antiparticles are less affected and their yield 
could also be well described with k=K2=2.0 GcV/fm 
(see also Fig. |3^,b). One note that very similar en- 
hancements are predicted for central and peripheral 
collisions, with shghtly higher value predicted for pe- 
ripheral (60-90%) coUisions, that could be attributed 
to different strength of quenching effect. 

In contrast to i?AuAu, ^cp obtained as scaled ratios 
of central (0-5%) and peripheral (60-90 %) Au+Au 
collisions show for (multi)strange particles a slight 
suppression relative to binary coUisions scaling (i?cp), 
see Fig. 3b, d in Ref. @. This suppression is mainly 
due to a final state effect, "jet quenching", which is 
stronger in central than in peripheral collisions. A 
striking difference between i?AuAu and i?cp has been 
reported by STAR [2l|. This experimental fact could 
be explained in our model as a consequence of inter- 
play of initial and final state effects. Baryon junc- 
tion loops (J J) and SCF effects are taken into consid- 
eration in both central and peripheral collisions. In 
p+p collisions the contribution of junction loops (J J) 
is negligible (due to small probability and kinemati- 
cal cuts), and SCF have a reduced strength (ki=1.5 
GeV/fm). This show that there could be a significant 
difference in the value and the meaning of i?cp and 
-RauAu due to differences in the baseline used for com- 
parison with binary scaling, i.e., peripheral (Au+Au) 
yields for Rqp and p+p yields for i?AuAu- 



IV. SUMMARY AND CONCLUSIONS 

We studied the influence of possible multi-gluon dy- 
namics (gluon junctions) and strong longitudinal color 
flclds (SCF) on particle production in heavy-ion colli- 
sions. A new parameterization of the p+p interaction 
based on new constituent and current quark masses 
for diquark is introduced. It leads to a simultaneously 
good description of non-strange and strange sector in 
p+p collisions. We show that J J loops play an im- 
portant role in particle production at mid-rapidity in 
d+Au and Au+Au colisions at RHIC energies. Intro- 
ducing a new J J loops algorithm in the framework of 
HIJING/BB v2.0, leads to a consistent and signiflcant 
improvement in the description of the recent experi- 
mental results for protons, pions, and kaons for both 
reactions. The present studies within our model are 
limited to the effect of initial state baryon produc- 
tion via possible junction dynamics in strong flelds. 
It would be very interesting to consider a generaliza- 
tion of back reaction effects [5(| to the case not only of 
pair production relevant for mesons but to the more 
difficult three string junction configurations needed to 
describe baryon production. Baryon productions via 
the conventional default quark-diquark mechanisms in 
the Lund string model are known to be inadequate 
even in e+ + e^ phenomenology. This is one of the 
main reasons for our continued exploration of baryon 
junction alternative mechanisms. 

The strange and multi-strange particles could only 
be described in the framework of string models if we 
consider SCF effects. The mechanisms of their pro- 
duction is very sensitive to the early phase of nu- 
clear collisions, when fiuctuations in the color field 
strength are highest. Within HIJING/BB v2.0, SCF 
effects are modeled by varying the effective string 
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tension that controls the qq and qqqq pair creation 
rates and strangeness suppression factors. The mid- 
rapidity yield of (niulti)strange particles favor a large 
value of the average string tension for both colli- 
sions (i.e., K « 2.0 GeV/fm in rf+Au, and k « 
3.0 GeV/fm in Au+Au). A strong enhancement 
in strange baryon nuclear modification factors Raa 
with increasing strangeness content is predicted for 
d-l-Au and Au-|-Au. In contrast, a clear ordering with 
strangeness content is not predicted for lower mean 
values of string tension. 

A greater sensitivity to SCF effects is predicted for 
the nuclear modification factors of multi-strange hy- 
perons S and O. In particular, the measurement of fl 
and Q yields would provide an important test of the 
relevance of SCF fluctuations, helping us to choose ap- 
propiate values for the suppression factors 7q (where 
Q=qq, us, s, u), which have strong dependence on 
the main parameters of QCD (the constituent and 
current quark masses) and on the system size. Even 
though the success of this procedure has been clearly 
illustrated here, the full understanding of the produc- 
tion of (multi)strange particles in relativistic heavy- 
ion collisions remain an exciting open question, and 
challenge many theoretical ideas. 
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APPENDIX: TABLES 



TABLE I: Main parameters used in the calculation for p+p 
collisions. The parameters are defined in the text. Set 1 is 
used in PYTHIA (and regular HIJING). Set 2 is obtained 
using the constituent quark masses from Sec. II B. Set 3 
includes an additional increase of the string tension to ki 
= 1.5 GeV/fm. 



p + p 


At (GeV/fm) 


7qq 


7s 


7us 


7io 


cr, (GeV/c) 


/ 


Set 1 


KO = 1.0 


0.09 


0.30 


0.40 


0.05 


0.350 


1 


Set 2 


KO = 1.0 


0.02 


0.30 


0.40 


0.05 


0.350 


1 


Set 3 


Kl = 1.5 


0.07 


0.45 


0.54 


0.09 


0.430 


1 



TABLE II: Main parameters used in d+Au and Au+Au 
collisions. The parameters are defined in the text. Set 
1 correspond to calculations without J J loops and SCF 
effects. Set 2 adds the contribution of J J loops. Set 3-5 
include both effects and correspond to different values of 
the string tension. 

(d)Au + Au K (GcV/fm) 7qq 7a 7us 710 aq (GeV/c) / 



Set 1 


KO 


= 1.0 


0.02 


0.30 


0.40 


0.05 


0.350 


Set 2 


KO 


= 1.0 


0.02 


0.30 


0.40 


0.05 


0.350 


Set 3 


Kl 


= 1.5 


0.07 


0.45 


0.54 


0.09 


0.430 


Set 4 


K2 


= 2.0 


0.14 


0.55 


0.63 


0.12 


0.495 


Set 5 




= 3.0 


0.27 


0.67 


0.74 


0.18 


0.606 
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